Analises 2022
Rafael Catoia
12/10/2022
library(dplyr) ; library(ggplot2) ; library(ggthemes) ; library(tidyr)## Functions
kable_cust = function(df_aux,ndec=5,cap=''){
return(df_aux %>% knitr::kable(digits = ndec,caption = cap) %>% kableExtra::kable_styling())
}Estudo
Novo relatório com análises para Ana Júlia - Los Hermanos
IRGA
UFSCar
CCSM
Rhizotron
dframe_rhiz = data.table::fread("/Users/rafaelcatoia/MyDrive/08_Consultorias/01_Ana_Julia/04_Projeto4/Dados/Rizhotron.txt")
dframe_rhiz_ufscar = dframe_rhiz %>%
filter(Intituto=='UFSCar') %>% select(-Manejo,-Tratamentos) %>%
mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')),
Rhizotron_2022=10*Rhizotron_2022) %>%
pivot_longer(cols = c("Rhizotron_2021","Rhizotron_2022"),names_to = 'Safra',
values_to = 'Rhizotron') %>%
mutate(Safra=gsub('Rhizotron_','',Safra))
dframe_rhiz_ccsm = dframe_rhiz %>%
filter(Intituto=='CCSM') %>% select(-Tratamento,-Tratamentos) %>%
mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV'))) %>%
pivot_longer(cols = c("Rhizotron_2021","Rhizotron_2022"),names_to = 'Safra',
values_to = 'Rhizotron') %>%
mutate(Safra=gsub('Rhizotron_','',Safra))UFSCar
dframe_summ_ufscar = dframe_rhiz_ufscar %>%
select(Safra,RegHidrico,Tratamento,Rhizotron) %>%
group_by(Safra,RegHidrico,Tratamento) %>%
summarise(Mean = mean(Rhizotron,na.rm=T),
SD = sd(Rhizotron,na.rm=T)) %>%
mutate(`Mean ± SD`=paste(round(Mean,5),'±',round(SD,5)))Graficos
Gráfico De Barras
Option 1
yname = 'Média de rhizotron'
dframe_summ_ufscar %>% ggplot(
aes(x=Tratamento,
fill=RegHidrico,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
facet_wrap(~Safra)+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(yname)
Option 2
dframe_summ_ufscar %>% ggplot(
aes(x=RegHidrico,
fill=Tratamento,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
facet_wrap(~Safra)+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(yname)
Gráfico De Barras com erro
Option 1
dframe_summ_ufscar %>% ggplot(
aes(x=Tratamento,
fill=RegHidrico,
y=Mean)) +
geom_bar(stat = 'identity',position = position_dodge(),alpha=0.9)+
geom_errorbar(
aes(x=Tratamento,y=Mean,
color=RegHidrico,
#linetype=RegHidrico,
ymin=Mean-SD,ymax=Mean+SD),
width = 0.2,alpha=0.8,stat = 'identity',color='black',
position = position_dodge(0.9))+
facet_wrap(~Safra)+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(yname)
Option 2
dframe_summ_ufscar %>% ggplot(
aes(x=RegHidrico,
fill=Tratamento,
y=Mean)) +
geom_bar(stat = 'identity',position = position_dodge(),alpha=0.9)+
geom_errorbar(
aes(x=RegHidrico,y=Mean,
color=Tratamento,
#linetype=RegHidrico,
ymin=Mean-SD,ymax=Mean+SD),
width = 0.2,alpha=0.8,stat = 'identity',color='black',
position = position_dodge(0.9))+
facet_wrap(~Safra)+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(yname)
Gráfico de Pontos + Média e Erro
Option 1
dframe_rhiz_ufscar %>% ggplot()+
facet_wrap(~Safra,scales='free')+
geom_point(aes(x=Tratamento,y=Rhizotron,
color=RegHidrico,
shape=RegHidrico),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(x=Tratamento,y=Mean,
color=RegHidrico,
shape=RegHidrico),
data=dframe_summ_ufscar,
position = position_dodge(width = .75))+
geom_errorbar(aes(x=Tratamento,y=Mean,
color=RegHidrico,
linetype=RegHidrico,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ_ufscar,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()
Option 2
dframe_rhiz_ufscar %>% ggplot()+
facet_wrap(~Safra,scales='free')+
geom_point(aes(x=RegHidrico,y=Rhizotron,
color=Tratamento,
shape=Tratamento),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(x=RegHidrico,y=Mean,
color=Tratamento,
shape=Tratamento),
data=dframe_summ_ufscar,
position = position_dodge(width = .75))+
geom_errorbar(aes(x=RegHidrico,y=Mean,
color=Tratamento,
linetype=Tratamento,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ_ufscar,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()
Tabela Resumo
dframe_summ_ufscar %>% select(Safra,RegHidrico,Tratamento,`Mean ± SD`) %>% kable_cust(cap = 'Média e Desvio Padrao - Rhizotron')| Safra | RegHidrico | Tratamento | Mean ± SD |
|---|---|---|---|
| 2021 | S | BRF_PF+FD | 0.01001 ± 0.00631 |
| 2021 | S | H10+FD | 0.0105 ± 0.00623 |
| 2021 | S | H10+H152 | 0.03834 ± 0.01675 |
| 2021 | TRF | BRF_PF+FD | 0.02045 ± 0.02714 |
| 2021 | TRF | H10+FD | 0.00918 ± 0.01088 |
| 2021 | TRF | H10+H152 | 0.01859 ± 0.0209 |
| 2021 | TRV | BRF_PF+FD | 0.00808 ± 0.00513 |
| 2021 | TRV | H10+FD | 0.01402 ± 0.0136 |
| 2021 | TRV | H10+H152 | 0.02275 ± 0.00243 |
| 2022 | S | BRF_PF+FD | 0.03457 ± 0.00308 |
| 2022 | S | H10+FD | 0.03134 ± 0.00449 |
| 2022 | S | H10+H152 | 0.04866 ± 0.00287 |
| 2022 | TRF | BRF_PF+FD | 0.02024 ± 0.00499 |
| 2022 | TRF | H10+FD | 0.01622 ± 0.00192 |
| 2022 | TRF | H10+H152 | 0.02094 ± 0.00723 |
| 2022 | TRV | BRF_PF+FD | 0.01795 ± 0.00295 |
| 2022 | TRV | H10+FD | 0.02425 ± 0.00266 |
| 2022 | TRV | H10+H152 | 0.02961 ± 0.00519 |
Modelo
2021
df_tk <- dframe_rhiz_ufscar %>% filter(Safra=='2021')
fit_lm <- lm(data=df_tk,
formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)
fit_aov <- aov(data=df_tk,
formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)
summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)
summ_tk <- df_tk %>% group_by(RegHidrico,Tratamento) %>%
summarise(mean=mean(Rhizotron)) %>% arrange(-mean)
summ_tk$Group =as.data.frame.list(mtcomp$`RegHidrico:Tratamento`)[,1]
summ_tk %>% kable_cust()| RegHidrico | Tratamento | mean | Group |
|---|---|---|---|
| S | H10+H152 | 0.03834 | a |
| TRV | H10+H152 | 0.02275 | a |
| TRF | BRF_PF+FD | 0.02045 | a |
| TRF | H10+H152 | 0.01859 | a |
| TRV | H10+FD | 0.01402 | a |
| S | H10+FD | 0.01050 | a |
| S | BRF_PF+FD | 0.01001 | a |
| TRF | H10+FD | 0.00918 | a |
| TRV | BRF_PF+FD | 0.00808 | a |
#tx <- with(data = df_tk,interaction(RegHidrico,Tratamento))
#amod <- aov(Rhizotron ~ RegHidrico+Tratamento+tx, data=df_tk)
#Tukey = agricolae::HSD.test(amod,trt = 'tx',group=TRUE)
#Tukey$groups %>% kable_cust()par(mfrow=c(2,2))
plot(fit_lm)
par(mfrow=c(1,1))summ_aov## Df Sum Sq Mean Sq F value Pr(>F)
## RegHidrico 2 0.000107 0.0000534 0.257 0.7761
## Tratamento 2 0.001277 0.0006387 3.077 0.0709 .
## RegHidrico:Tratamento 4 0.000847 0.0002116 1.020 0.4236
## Residuals 18 0.003736 0.0002076
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1summ_lm$coefficients %>% kable_cust()| Estimate | Std. Error | t value | Pr(>|t|) | |
|---|---|---|---|---|
| (Intercept) | 0.01001 | 0.00832 | 1.20307 | 0.24454 |
| RegHidricoTRF | 0.01044 | 0.01176 | 0.88738 | 0.38658 |
| RegHidricoTRV | -0.00192 | 0.01176 | -0.16355 | 0.87191 |
| TratamentoH10+FD | 0.00049 | 0.01176 | 0.04153 | 0.96733 |
| TratamentoH10+H152 | 0.02833 | 0.01176 | 2.40859 | 0.02695 |
| RegHidricoTRF:TratamentoH10+FD | -0.01176 | 0.01664 | -0.70673 | 0.48878 |
| RegHidricoTRV:TratamentoH10+FD | 0.00544 | 0.01664 | 0.32723 | 0.74727 |
| RegHidricoTRF:TratamentoH10+H152 | -0.03019 | 0.01664 | -1.81452 | 0.08630 |
| RegHidricoTRV:TratamentoH10+H152 | -0.01367 | 0.01664 | -0.82149 | 0.42211 |
2022
df_tk <- dframe_rhiz_ufscar %>% filter(Safra=='2022')
fit_lm <- lm(data=df_tk,
formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)
fit_aov <- aov(data=df_tk,
formula = Rhizotron~RegHidrico+Tratamento+RegHidrico*Tratamento)
summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)
summ_tk <- df_tk %>% group_by(RegHidrico,Tratamento) %>%
summarise(mean=mean(Rhizotron)) %>% arrange(-mean)
summ_tk$Group = as.data.frame.list(mtcomp$`RegHidrico:Tratamento`)[,1]
summ_tk %>% kable_cust()| RegHidrico | Tratamento | mean | Group |
|---|---|---|---|
| S | H10+H152 | 0.04866 | a |
| S | BRF_PF+FD | 0.03457 | b |
| S | H10+FD | 0.03134 | bc |
| TRV | H10+H152 | 0.02961 | bcd |
| TRV | H10+FD | 0.02425 | bcde |
| TRF | H10+H152 | 0.02094 | cde |
| TRF | BRF_PF+FD | 0.02024 | cde |
| TRV | BRF_PF+FD | 0.01795 | de |
| TRF | H10+FD | 0.01622 | e |
#tx <- with(data = df_tk,interaction(RegHidrico,Tratamento))
#amod <- aov(Rhizotron ~ RegHidrico+Tratamento+tx, data=df_tk)
#
#Tukey = agricolae::HSD.test(amod,trt = 'tx',group=TRUE)
#Tukey$groups %>% kable_cust()par(mfrow=c(2,2))
plot(fit_lm)
par(mfrow=c(1,1))summ_aov## Df Sum Sq Mean Sq F value Pr(>F)
## RegHidrico 2 0.0017678 0.0008839 49.258 5.01e-08 ***
## Tratamento 2 0.0004839 0.0002419 13.483 0.000264 ***
## RegHidrico:Tratamento 4 0.0002684 0.0000671 3.739 0.021999 *
## Residuals 18 0.0003230 0.0000179
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1summ_lm$coefficients %>% kable_cust()| Estimate | Std. Error | t value | Pr(>|t|) | |
|---|---|---|---|---|
| (Intercept) | 0.03457 | 0.00245 | 14.13354 | 0.00000 |
| RegHidricoTRF | -0.01433 | 0.00346 | -4.14327 | 0.00061 |
| RegHidricoTRV | -0.01661 | 0.00346 | -4.80345 | 0.00014 |
| TratamentoH10+FD | -0.00323 | 0.00346 | -0.93337 | 0.36298 |
| TratamentoH10+H152 | 0.01409 | 0.00346 | 4.07497 | 0.00071 |
| RegHidricoTRF:TratamentoH10+FD | -0.00079 | 0.00489 | -0.16097 | 0.87391 |
| RegHidricoTRV:TratamentoH10+FD | 0.00953 | 0.00489 | 1.94779 | 0.06721 |
| RegHidricoTRF:TratamentoH10+H152 | -0.01339 | 0.00489 | -2.73656 | 0.01356 |
| RegHidricoTRV:TratamentoH10+H152 | -0.00244 | 0.00489 | -0.49902 | 0.62381 |
CCSM
dframe_summ_ccsm = dframe_rhiz_ccsm %>%
select(Safra,Manejo,Copa,PortaEnxerto,Rhizotron) %>%
group_by(Safra,Manejo,Copa,PortaEnxerto) %>%
summarise(Mean = mean(Rhizotron,na.rm=T),
SD = sd(Rhizotron,na.rm=T)) %>%
mutate(`Mean ± SD`=paste(round(Mean,5),'±',round(SD,5)))Graficos
Gráfico de Barras
Option 1
yname = 'Média de rhizotron'
dframe_summ_ccsm %>% ggplot(
aes(x=PortaEnxerto,
fill=Copa,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
facet_wrap(Safra~Manejo,nrow=1)+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(yname)
Option 2
dframe_summ_ccsm %>% ggplot(
aes(x=Copa,
fill=PortaEnxerto,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
facet_wrap(Safra~Manejo,nrow=1)+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(yname)
Gráfico de Barras com erro
Option 1
dframe_summ_ccsm %>% ggplot(
aes(x=PortaEnxerto,
fill=Copa,
y=Mean)) +
geom_bar(stat = 'identity',position = position_dodge(),alpha=0.9)+
geom_errorbar(
aes(x=PortaEnxerto,y=Mean,
color=Copa,
ymin=Mean-SD,ymax=Mean+SD),
width = 0.2,alpha=0.8,stat = 'identity',color='black',
position = position_dodge(0.9))+
facet_wrap(Safra~Manejo,nrow=1)+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(yname)
Option 2
dframe_summ_ccsm %>% ggplot(
aes(x=Copa,
fill=PortaEnxerto,
y=Mean)) +
geom_bar(stat = 'identity',position = position_dodge(),alpha=0.9)+
geom_errorbar(
aes(x=Copa,y=Mean,
color=PortaEnxerto,
ymin=Mean-SD,ymax=Mean+SD),
width = 0.2,alpha=0.8,stat = 'identity',color='black',
position = position_dodge(0.9))+
facet_wrap(Safra~Manejo,nrow=1)+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(yname)
Gráfico de Pontos + Média e Erro
Option 1
dframe_rhiz_ccsm %>% ggplot()+
facet_wrap(~Safra)+
geom_point(aes(x=PortaEnxerto,
y=Rhizotron,
color=Copa,
shape=Manejo),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(x=PortaEnxerto,
y=Mean,
color=Copa,
shape=Manejo),
data=dframe_summ_ccsm,
position = position_dodge(width = .75))+
geom_errorbar(aes(x=PortaEnxerto,
y=Mean,
color=Copa,
shape=Manejo,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ_ccsm,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()
Option 2
dframe_rhiz_ccsm %>% ggplot()+
facet_wrap(~Safra)+
geom_point(aes(x=Copa,
y=Rhizotron,
color=PortaEnxerto,
shape=Manejo),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(x=Copa,
y=Mean,
color=PortaEnxerto,
shape=Manejo),
data=dframe_summ_ccsm,
position = position_dodge(width = .75))+
geom_errorbar(aes(x=Copa,
y=Mean,
color=PortaEnxerto,
shape=Manejo,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ_ccsm,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()
Tabela Resumo
dframe_summ_ccsm %>% select(Safra,Manejo,PortaEnxerto,Copa,`Mean ± SD`) %>% kable_cust(cap = 'Média e Desvio Padrao - Rhizotron')| Safra | Manejo | PortaEnxerto | Copa | Mean ± SD |
|---|---|---|---|---|
| 2021 | Convencional | FD | BRF_PF | 0.00521 ± 0.00401 |
| 2021 | Convencional | H152 | BRF_PF | 0.00837 ± 0.0044 |
| 2021 | Convencional | FD | IAC10 | 0.0035 ± 0.00158 |
| 2021 | Convencional | H152 | IAC10 | 0.01423 ± 0.00303 |
| 2021 | Convencional | FD | IAC5 | 0.01404 ± 0.01049 |
| 2021 | Convencional | H152 | IAC5 | 0.00987 ± 0.00401 |
| 2021 | Ecologico | FD | BRF_PF | 0.0088 ± 0.0014 |
| 2021 | Ecologico | H152 | BRF_PF | 0.01557 ± 0.00379 |
| 2021 | Ecologico | FD | IAC10 | 0.00549 ± 0.00033 |
| 2021 | Ecologico | H152 | IAC10 | 0.03573 ± 0.01341 |
| 2021 | Ecologico | FD | IAC5 | 0.01022 ± 0.01135 |
| 2021 | Ecologico | H152 | IAC5 | 0.01021 ± 0.00289 |
| 2022 | Convencional | FD | BRF_PF | 0.01011 ± 0.00569 |
| 2022 | Convencional | H152 | BRF_PF | 0.02145 ± 0.00038 |
| 2022 | Convencional | FD | IAC10 | 0.00386 ± 0.00174 |
| 2022 | Convencional | H152 | IAC10 | 0.00909 ± 0.00083 |
| 2022 | Convencional | FD | IAC5 | 0.00171 ± 6e-05 |
| 2022 | Convencional | H152 | IAC5 | 0.01492 ± 0.00074 |
| 2022 | Ecologico | FD | BRF_PF | 0.00425 ± 0.00071 |
| 2022 | Ecologico | H152 | BRF_PF | 0.02376 ± 0.00298 |
| 2022 | Ecologico | FD | IAC10 | 0.00566 ± 0.00022 |
| 2022 | Ecologico | H152 | IAC10 | 0.00891 ± 0.00142 |
| 2022 | Ecologico | FD | IAC5 | 0.00459 ± 0.00013 |
| 2022 | Ecologico | H152 | IAC5 | 0.01106 ± 0.00222 |
Modelo
2021
df_tk <- dframe_rhiz_ccsm %>% filter(Safra=='2021')
fit_lm <- lm(data=df_tk,
formula = Rhizotron~Manejo+Copa +PortaEnxerto+
Manejo*Copa+Manejo*PortaEnxerto +
Copa*PortaEnxerto+
Manejo*Copa*PortaEnxerto
)
fit_aov <- aov(data=df_tk,
formula = Rhizotron~Manejo+Copa +PortaEnxerto+
Manejo*Copa+#Manejo*PortaEnxerto +
Copa*PortaEnxerto+#+
Manejo*Copa*PortaEnxerto
)
summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)
summ_tk <- df_tk %>% group_by(Manejo,Copa,PortaEnxerto) %>%
summarise(mean=mean(Rhizotron)) %>% arrange(-mean)
summ_tk$Group =as.data.frame.list(mtcomp$`Manejo:Copa:PortaEnxerto`)[,1]
summ_tk %>% kable_cust()| Manejo | Copa | PortaEnxerto | mean | Group |
|---|---|---|---|---|
| Ecologico | IAC10 | H152 | 0.03573 | a |
| Ecologico | BRF_PF | H152 | 0.01557 | b |
| Convencional | IAC10 | H152 | 0.01423 | b |
| Convencional | IAC5 | FD | 0.01404 | b |
| Ecologico | IAC5 | FD | 0.01022 | b |
| Ecologico | IAC5 | H152 | 0.01021 | b |
| Convencional | IAC5 | H152 | 0.00987 | b |
| Ecologico | BRF_PF | FD | 0.00880 | b |
| Convencional | BRF_PF | H152 | 0.00837 | b |
| Ecologico | IAC10 | FD | 0.00549 | b |
| Convencional | BRF_PF | FD | 0.00521 | b |
| Convencional | IAC10 | FD | 0.00350 | b |
par(mfrow=c(2,2))
plot(fit_lm)
par(mfrow=c(1,1))summ_aov## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0002372 0.0002372 5.620 0.026137 *
## Copa 2 0.0001736 0.0000868 2.056 0.149916
## PortaEnxerto 1 0.0005457 0.0005457 12.928 0.001454 **
## Manejo:Copa 2 0.0002731 0.0001366 3.236 0.057000 .
## Copa:PortaEnxerto 2 0.0007999 0.0004000 9.476 0.000926 ***
## Manejo:PortaEnxerto 1 0.0001861 0.0001861 4.409 0.046447 *
## Manejo:Copa:PortaEnxerto 2 0.0001221 0.0000611 1.447 0.255160
## Residuals 24 0.0010130 0.0000422
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1summ_lm$coefficients %>% kable_cust()| Estimate | Std. Error | t value | Pr(>|t|) | |
|---|---|---|---|---|
| (Intercept) | 0.00521 | 0.00375 | 1.38946 | 0.17745 |
| ManejoEcologico | 0.00359 | 0.00530 | 0.67625 | 0.50535 |
| CopaIAC10 | -0.00171 | 0.00530 | -0.32305 | 0.74946 |
| CopaIAC5 | 0.00883 | 0.00530 | 1.66454 | 0.10901 |
| PortaEnxertoH152 | 0.00316 | 0.00530 | 0.59530 | 0.55722 |
| ManejoEcologico:CopaIAC10 | -0.00160 | 0.00750 | -0.21266 | 0.83339 |
| ManejoEcologico:CopaIAC5 | -0.00741 | 0.00750 | -0.98735 | 0.33332 |
| ManejoEcologico:PortaEnxertoH152 | 0.00362 | 0.00750 | 0.48201 | 0.63416 |
| CopaIAC10:PortaEnxertoH152 | 0.00757 | 0.00750 | 1.00905 | 0.32302 |
| CopaIAC5:PortaEnxertoH152 | -0.00733 | 0.00750 | -0.97650 | 0.33856 |
| ManejoEcologico:CopaIAC10:PortaEnxertoH152 | 0.01589 | 0.01061 | 1.49799 | 0.14717 |
| ManejoEcologico:CopaIAC5:PortaEnxertoH152 | 0.00054 | 0.01061 | 0.05116 | 0.95962 |
2022
df_tk <- dframe_rhiz_ccsm %>% filter(Safra=='2022')
fit_lm <- lm(data=df_tk,
formula = Rhizotron~
#Manejo+
Copa +PortaEnxerto+
#Manejo*Copa+Manejo*PortaEnxerto +
Copa*PortaEnxerto#+
#Manejo*Copa*PortaEnxerto
)
fit_aov <- aov(data=df_tk,
formula = Rhizotron~
Manejo+Copa +PortaEnxerto+
Manejo*Copa+Manejo*PortaEnxerto +
Copa*PortaEnxerto+
Manejo*Copa*PortaEnxerto
)
summ_lm = summary(fit_lm)
summ_aov = summary(fit_aov)tk <- TukeyHSD(fit_aov)
mtcomp <- multcompView::multcompLetters4(fit_aov,tk)
summ_tk <- df_tk %>% group_by(Copa,PortaEnxerto) %>%
summarise(mean=mean(Rhizotron)) %>% arrange(-mean)
summ_tk$Group =as.data.frame.list(mtcomp$`Copa:PortaEnxerto`)[,1]
summ_tk %>% kable_cust()| Copa | PortaEnxerto | mean | Group |
|---|---|---|---|
| BRF_PF | H152 | 0.02261 | a |
| IAC5 | H152 | 0.01299 | b |
| IAC10 | H152 | 0.00900 | c |
| BRF_PF | FD | 0.00718 | cd |
| IAC10 | FD | 0.00476 | de |
| IAC5 | FD | 0.00315 | e |
par(mfrow=c(2,2))
plot(fit_lm)
par(mfrow=c(1,1))
par(mfrow=c(2,2))
plot(fit_aov)
par(mfrow=c(1,1))summ_aov## Df Sum Sq Mean Sq F value Pr(>F)
## Manejo 1 0.0000021 0.0000021 0.480 0.495002
## Copa 2 0.0004488 0.0002244 50.641 2.44e-09 ***
## PortaEnxerto 1 0.0008701 0.0008701 196.364 4.73e-13 ***
## Manejo:Copa 2 0.0000099 0.0000050 1.122 0.342089
## Manejo:PortaEnxerto 1 0.0000001 0.0000001 0.017 0.896199
## Copa:PortaEnxerto 2 0.0001877 0.0000938 21.174 5.02e-06 ***
## Manejo:Copa:PortaEnxerto 2 0.0000871 0.0000435 9.827 0.000763 ***
## Residuals 24 0.0001063 0.0000044
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1summ_lm$coefficients %>% kable_cust()| Estimate | Std. Error | t value | Pr(>|t|) | |
|---|---|---|---|---|
| (Intercept) | 0.00718 | 0.00107 | 6.71937 | 0.00000 |
| CopaIAC10 | -0.00242 | 0.00151 | -1.60200 | 0.11963 |
| CopaIAC5 | -0.00403 | 0.00151 | -2.66480 | 0.01228 |
| PortaEnxertoH152 | 0.01543 | 0.00151 | 10.20594 | 0.00000 |
| CopaIAC10:PortaEnxertoH152 | -0.01118 | 0.00214 | -5.23293 | 0.00001 |
| CopaIAC5:PortaEnxertoH152 | -0.00559 | 0.00214 | -2.61646 | 0.01378 |
Potencial Hidrico
UFSCar
CCSM
Desenvolvimento Vegetativo
dframe_dv = data.table::fread("/Users/rafaelcatoia/MyDrive/08_Consultorias/01_Ana_Julia/04_Projeto4/Dados/DesenvVegetativo.txt")
#dframe_qf_ufscar = dframe_qf %>%
# filter(Instituto=='UFSCar') %>% select(-Manejo) %>%
# mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))
#
#dframe_qf_ccsm = dframe_qf %>%
# filter(Instituto=='CCSM') %>% select(-Tratamento) %>%
# mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))UFSCar
CCSM
Qualidade do fruto
dframe_qf = data.table::fread("/Users/rafaelcatoia/MyDrive/08_Consultorias/01_Ana_Julia/04_Projeto4/Dados/QualidadeFruto.txt")
dframe_qf_ufscar = dframe_qf %>%
filter(Instituto=='UFSCar') %>% select(-Manejo) %>%
mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))
dframe_qf_ccsm = dframe_qf %>%
filter(Instituto=='CCSM') %>% select(-Tratamento) %>%
mutate(RegHidrico =factor(RegHidrico,levels = c('S','TRF','TRV')))summ_qf <- function(dframe,variable_name = ''){
if (variable_name==''){
stop()
}
dframe = dframe %>% select(one_of(variable_name),RegHidrico,Tratamento)
colnames(dframe)[1]<-'Y'
out<- list()
dframe_summ <- dframe %>% group_by(RegHidrico,Tratamento) %>%
summarise(Mean = mean(Y,na.rm=T),
SD = sd(Y,na.rm=T)) %>%
mutate(`Mean ± SD`=paste(round(Mean,2),'±',round(SD,2)))
out$dframe_summ <- dframe_summ
out$barp <- dframe_summ %>% ggplot(
aes(x=RegHidrico,
fill=Tratamento,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$barp_2 <- dframe_summ %>% ggplot(
aes(x=Tratamento,
fill=RegHidrico,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$barp_error <- dframe_summ %>% ggplot(
aes(x=RegHidrico,
fill=Tratamento,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
geom_errorbar(
aes(x=RegHidrico,y=Mean,
color=Tratamento,
ymin=Mean-SD,ymax=Mean+SD),
width = 0.2,alpha=0.8,stat = 'identity',color='black',
position = position_dodge(0.9))+
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$barp_error_2 <- dframe_summ %>% ggplot(
aes(x=Tratamento,
fill=RegHidrico,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
geom_errorbar(
aes(x=Tratamento,y=Mean,
color=RegHidrico,
ymin=Mean-SD,ymax=Mean+SD),
width = 0.2,alpha=0.8,stat = 'identity',color='black',
position = position_dodge(0.9))+
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$mean_desv_plot <-
dframe %>% ggplot()+
geom_point(aes(
x=RegHidrico,
y=Y,
color=Tratamento,
shape=Tratamento),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(
x=RegHidrico,
y=Mean,
color=Tratamento,
shape=Tratamento),
data=dframe_summ,
position = position_dodge(width = .75))+
geom_errorbar(aes(
x=RegHidrico,
y=Mean,
color=Tratamento,
shape=Tratamento,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$mean_desv_plot_2 <-
dframe %>% ggplot()+
geom_point(aes(
x=Tratamento,
y=Y,
color=RegHidrico,
shape=RegHidrico),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(
x=Tratamento,
y=Mean,
color=RegHidrico,
shape=RegHidrico),
data=dframe_summ,
position = position_dodge(width = .75))+
geom_errorbar(aes(
x=Tratamento,
y=Mean,
color=RegHidrico,
shape=RegHidrico,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
return(out)
}
summ_qf_ccsm <- function(dframe,variable_name = ''){
if (variable_name==''){
stop()
}
dframe = dframe %>% select(one_of(variable_name),Manejo,Copa,PortaEnxerto)
colnames(dframe)[1]<-'Y'
out<- list()
dframe = dframe %>% mutate(Tratamento = paste(Copa,'+',PortaEnxerto,sep=''))
dframe_summ <- dframe %>% group_by(Manejo,Tratamento) %>%
summarise(Mean = mean(Y,na.rm=T),
SD = sd(Y,na.rm=T)) %>%
mutate(`Mean ± SD`=paste(round(Mean,2),'±',round(SD,2)))
out$dframe_summ <- dframe_summ
out$barp <- dframe_summ %>% ggplot(
aes(x=Manejo,
fill=Tratamento,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$barp_2 <- dframe_summ %>% ggplot(
aes(x=Tratamento,
fill=Manejo,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$barp_error <- dframe_summ %>% ggplot(
aes(x=Manejo,
fill=Tratamento,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
geom_errorbar(
aes(x=Manejo,y=Mean,
color=Tratamento,
ymin=Mean-SD,ymax=Mean+SD),
width = 0.2,alpha=0.8,stat = 'identity',color='black',
position = position_dodge(0.9))+
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$barp_error_2 <- dframe_summ %>% ggplot(
aes(x=Tratamento,
fill=Manejo,
y=Mean)) +
geom_bar(stat = 'identity',position = 'dodge')+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
geom_errorbar(
aes(x=Tratamento,y=Mean,
color=Manejo,
ymin=Mean-SD,ymax=Mean+SD),
width = 0.2,alpha=0.8,stat = 'identity',color='black',
position = position_dodge(0.9))+
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$mean_desv_plot <-
dframe %>% ggplot()+
geom_point(aes(
x=Manejo,
y=Y,
color=Tratamento,
shape=Tratamento),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(
x=Manejo,
y=Mean,
color=Tratamento,
shape=Tratamento),
data=dframe_summ,
position = position_dodge(width = .75))+
geom_errorbar(aes(
x=Manejo,
y=Mean,
color=Tratamento,
shape=Tratamento,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
out$mean_desv_plot_2 <-
dframe %>% ggplot()+
geom_point(aes(
x=Tratamento,
y=Y,
color=Manejo,
shape=Manejo),
position = position_dodge(width = .75),alpha=0.5)+
geom_point(aes(
x=Tratamento,
y=Mean,
color=Manejo,
shape=Manejo),
data=dframe_summ,
position = position_dodge(width = .75))+
geom_errorbar(aes(
x=Tratamento,
y=Mean,
color=Manejo,
shape=Manejo,
ymin=Mean-SD,ymax=Mean+SD),
data=dframe_summ,
width = 0.2,alpha=0.9,
position = position_dodge(0.75))+
theme_minimal(base_size=14)+theme(legend.position = 'bottom') +
scale_color_gdocs()+
scale_fill_gdocs()+
ylab(paste('Média de', variable_name ))
return(out)
}UFScar
Massa
results <- summ_qf(dframe_qf_ufscar,variable_name = "Massa")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| RegHidrico | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| S | FD+BRS_PF | 81.66667 | 7.09460 | 81.67 ± 7.09 |
| S | FD+H10 | 47.00000 | 2.64575 | 47 ± 2.65 |
| S | H14+H10 | 69.33333 | 3.51188 | 69.33 ± 3.51 |
| S | H152+H10 | 69.33333 | 6.11010 | 69.33 ± 6.11 |
| S | H228+H10 | 47.66667 | 0.57735 | 47.67 ± 0.58 |
| TRF | FD+BRS_PF | 86.66667 | 8.14453 | 86.67 ± 8.14 |
| TRF | FD+H10 | 82.06667 | 4.10041 | 82.07 ± 4.1 |
| TRF | H14+H10 | 76.00000 | 2.64575 | 76 ± 2.65 |
| TRF | H152+H10 | 79.50000 | 8.50000 | 79.5 ± 8.5 |
| TRF | H228+H10 | 83.00000 | 8.18535 | 83 ± 8.19 |
| TRV | FD+BRS_PF | 71.66667 | 1.15470 | 71.67 ± 1.15 |
| TRV | FD+H10 | 75.33333 | 1.15470 | 75.33 ± 1.15 |
| TRV | H14+H10 | 54.66667 | 8.50490 | 54.67 ± 8.5 |
| TRV | H152+H10 | 74.66667 | 3.51188 | 74.67 ± 3.51 |
| TRV | H228+H10 | 68.33333 | 7.02377 | 68.33 ± 7.02 |
Altura
results <- summ_qf(dframe_qf_ufscar,variable_name = "Altura")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| RegHidrico | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| S | FD+BRS_PF | 5.80000 | 0.20000 | 5.8 ± 0.2 |
| S | FD+H10 | 4.86667 | 0.11547 | 4.87 ± 0.12 |
| S | H14+H10 | 5.53333 | 0.05774 | 5.53 ± 0.06 |
| S | H152+H10 | 5.40000 | 0.20000 | 5.4 ± 0.2 |
| S | H228+H10 | 5.06667 | 0.05774 | 5.07 ± 0.06 |
| TRF | FD+BRS_PF | 5.86667 | 0.28868 | 5.87 ± 0.29 |
| TRF | FD+H10 | 5.83333 | 0.15275 | 5.83 ± 0.15 |
| TRF | H14+H10 | 5.60000 | 0.00000 | 5.6 ± 0 |
| TRF | H152+H10 | 5.70000 | 0.10000 | 5.7 ± 0.1 |
| TRF | H228+H10 | 5.73333 | 0.25166 | 5.73 ± 0.25 |
| TRV | FD+BRS_PF | 5.43333 | 0.11547 | 5.43 ± 0.12 |
| TRV | FD+H10 | 5.50000 | 0.10000 | 5.5 ± 0.1 |
| TRV | H14+H10 | 5.13333 | 0.11547 | 5.13 ± 0.12 |
| TRV | H152+H10 | 5.46667 | 0.32146 | 5.47 ± 0.32 |
| TRV | H228+H10 | 5.53333 | 0.28868 | 5.53 ± 0.29 |
Largura
results <- summ_qf(dframe_qf_ufscar,variable_name = "Largura")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| RegHidrico | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| S | FD+BRS_PF | 5.26667 | 0.11547 | 5.27 ± 0.12 |
| S | FD+H10 | 4.36667 | 0.05774 | 4.37 ± 0.06 |
| S | H14+H10 | 5.00000 | 0.10000 | 5 ± 0.1 |
| S | H152+H10 | 5.03333 | 0.15275 | 5.03 ± 0.15 |
| S | H228+H10 | 4.40000 | 0.10000 | 4.4 ± 0.1 |
| TRF | FD+BRS_PF | 5.43333 | 0.15275 | 5.43 ± 0.15 |
| TRF | FD+H10 | 5.33333 | 0.11547 | 5.33 ± 0.12 |
| TRF | H14+H10 | 5.23333 | 0.05774 | 5.23 ± 0.06 |
| TRF | H152+H10 | 5.40000 | 0.10000 | 5.4 ± 0.1 |
| TRF | H228+H10 | 5.30000 | 0.17321 | 5.3 ± 0.17 |
| TRV | FD+BRS_PF | 5.10000 | 0.00000 | 5.1 ± 0 |
| TRV | FD+H10 | 5.16667 | 0.05774 | 5.17 ± 0.06 |
| TRV | H14+H10 | 4.63333 | 0.30551 | 4.63 ± 0.31 |
| TRV | H152+H10 | 4.93333 | 0.28868 | 4.93 ± 0.29 |
| TRV | H228+H10 | 4.96667 | 0.15275 | 4.97 ± 0.15 |
AL
results <- summ_qf(dframe_qf_ufscar,variable_name = "AL")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| RegHidrico | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| S | FD+BRS_PF | 1.10333 | 0.02082 | 1.1 ± 0.02 |
| S | FD+H10 | 1.11667 | 0.02517 | 1.12 ± 0.03 |
| S | H14+H10 | 1.10667 | 0.01155 | 1.11 ± 0.01 |
| S | H152+H10 | 1.07333 | 0.03055 | 1.07 ± 0.03 |
| S | H228+H10 | 1.15333 | 0.04041 | 1.15 ± 0.04 |
| TRF | FD+BRS_PF | 1.08333 | 0.02517 | 1.08 ± 0.03 |
| TRF | FD+H10 | 1.09333 | 0.02082 | 1.09 ± 0.02 |
| TRF | H14+H10 | 1.07333 | 0.01155 | 1.07 ± 0.01 |
| TRF | H152+H10 | 1.05667 | 0.00577 | 1.06 ± 0.01 |
| TRF | H228+H10 | 1.08333 | 0.02082 | 1.08 ± 0.02 |
| TRV | FD+BRS_PF | 1.06667 | 0.02309 | 1.07 ± 0.02 |
| TRV | FD+H10 | 1.06667 | 0.02309 | 1.07 ± 0.02 |
| TRV | H14+H10 | 1.11000 | 0.05000 | 1.11 ± 0.05 |
| TRV | H152+H10 | 1.11000 | 0.01000 | 1.11 ± 0.01 |
| TRV | H228+H10 | 1.11333 | 0.03055 | 1.11 ± 0.03 |
Suco
results <- summ_qf(dframe_qf_ufscar,variable_name = "Suco")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| RegHidrico | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| S | FD+BRS_PF | 48.76667 | 4.05010 | 48.77 ± 4.05 |
| S | FD+H10 | 37.46667 | 3.35012 | 37.47 ± 3.35 |
| S | H14+H10 | 45.63333 | 1.50444 | 45.63 ± 1.5 |
| S | H152+H10 | 45.13333 | 2.58908 | 45.13 ± 2.59 |
| S | H228+H10 | 39.86667 | 0.46188 | 39.87 ± 0.46 |
| TRF | FD+BRS_PF | 51.83333 | 2.10792 | 51.83 ± 2.11 |
| TRF | FD+H10 | 51.20000 | 2.07846 | 51.2 ± 2.08 |
| TRF | H14+H10 | 46.46667 | 1.76163 | 46.47 ± 1.76 |
| TRF | H152+H10 | 45.00000 | 2.70000 | 45 ± 2.7 |
| TRF | H228+H10 | 45.63333 | 3.13900 | 45.63 ± 3.14 |
| TRV | FD+BRS_PF | 45.13333 | 0.05774 | 45.13 ± 0.06 |
| TRV | FD+H10 | 49.56667 | 0.75056 | 49.57 ± 0.75 |
| TRV | H14+H10 | 38.63333 | 4.06735 | 38.63 ± 4.07 |
| TRV | H152+H10 | 41.90000 | 1.60935 | 41.9 ± 1.61 |
| TRV | H228+H10 | 43.43333 | 3.76475 | 43.43 ± 3.76 |
Acidez
results <- summ_qf(dframe_qf_ufscar,variable_name = "Acidez")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| RegHidrico | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| S | FD+BRS_PF | 5.29067 | 0.38493 | 5.29 ± 0.38 |
| S | FD+H10 | 5.42190 | 0.28282 | 5.42 ± 0.28 |
| S | H14+H10 | 5.01870 | 0.39804 | 5.02 ± 0.4 |
| S | H152+H10 | 4.60677 | 0.51469 | 4.61 ± 0.51 |
| S | H228+H10 | 5.06987 | 0.18971 | 5.07 ± 0.19 |
| TRF | FD+BRS_PF | 5.46163 | 0.26236 | 5.46 ± 0.26 |
| TRF | FD+H10 | 5.14260 | 0.43562 | 5.14 ± 0.44 |
| TRF | H14+H10 | 4.97363 | 0.08585 | 4.97 ± 0.09 |
| TRF | H152+H10 | 5.16717 | 0.15245 | 5.17 ± 0.15 |
| TRF | H228+H10 | 5.51613 | 0.34844 | 5.52 ± 0.35 |
| TRV | FD+BRS_PF | 5.15400 | 0.32922 | 5.15 ± 0.33 |
| TRV | FD+H10 | 4.99920 | 0.20471 | 5 ± 0.2 |
| TRV | H14+H10 | 5.00657 | 0.06020 | 5.01 ± 0.06 |
| TRV | H152+H10 | 4.79327 | 0.15489 | 4.79 ± 0.15 |
| TRV | H228+H10 | 5.50537 | 0.33448 | 5.51 ± 0.33 |
BrixCorr
results <- summ_qf(dframe_qf_ufscar,variable_name = "BrixCorr")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| RegHidrico | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| S | FD+BRS_PF | 8.90000 | 0.17321 | 8.9 ± 0.17 |
| S | FD+H10 | 8.80000 | 0.10000 | 8.8 ± 0.1 |
| S | H14+H10 | 8.83333 | 0.15275 | 8.83 ± 0.15 |
| S | H152+H10 | 8.63333 | 0.05774 | 8.63 ± 0.06 |
| S | H228+H10 | 9.56667 | 0.20817 | 9.57 ± 0.21 |
| TRF | FD+BRS_PF | 9.13333 | 0.15275 | 9.13 ± 0.15 |
| TRF | FD+H10 | 9.16667 | 0.05774 | 9.17 ± 0.06 |
| TRF | H14+H10 | 8.93333 | 0.11547 | 8.93 ± 0.12 |
| TRF | H152+H10 | 8.90000 | 0.10000 | 8.9 ± 0.1 |
| TRF | H228+H10 | 8.90000 | 0.00000 | 8.9 ± 0 |
| TRV | FD+BRS_PF | 8.83333 | 0.20817 | 8.83 ± 0.21 |
| TRV | FD+H10 | 9.20000 | 0.10000 | 9.2 ± 0.1 |
| TRV | H14+H10 | 8.86667 | 0.20817 | 8.87 ± 0.21 |
| TRV | H152+H10 | 8.66667 | 0.05774 | 8.67 ± 0.06 |
| TRV | H228+H10 | 8.96667 | 0.25166 | 8.97 ± 0.25 |
Ratio
results <- summ_qf(dframe_qf_ufscar,variable_name = "Ratio")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| RegHidrico | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| S | FD+BRS_PF | 1.70000 | 0.10000 | 1.7 ± 0.1 |
| S | FD+H10 | 1.60000 | 0.10000 | 1.6 ± 0.1 |
| S | H14+H10 | 1.73333 | 0.11547 | 1.73 ± 0.12 |
| S | H152+H10 | 1.90000 | 0.20000 | 1.9 ± 0.2 |
| S | H228+H10 | 1.90000 | 0.10000 | 1.9 ± 0.1 |
| TRF | FD+BRS_PF | 1.66667 | 0.05774 | 1.67 ± 0.06 |
| TRF | FD+H10 | 1.76667 | 0.15275 | 1.77 ± 0.15 |
| TRF | H14+H10 | 1.80000 | 0.00000 | 1.8 ± 0 |
| TRF | H152+H10 | 1.75000 | 0.05000 | 1.75 ± 0.05 |
| TRF | H228+H10 | 1.60000 | 0.10000 | 1.6 ± 0.1 |
| TRV | FD+BRS_PF | 1.73333 | 0.15275 | 1.73 ± 0.15 |
| TRV | FD+H10 | 1.86667 | 0.05774 | 1.87 ± 0.06 |
| TRV | H14+H10 | 1.76667 | 0.05774 | 1.77 ± 0.06 |
| TRV | H152+H10 | 1.80000 | 0.10000 | 1.8 ± 0.1 |
| TRV | H228+H10 | 1.63333 | 0.05774 | 1.63 ± 0.06 |
IT
results <- summ_qf(dframe_qf_ufscar,variable_name = "IT")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| RegHidrico | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| S | FD+BRS_PF | 1.77000 | 0.18358 | 1.77 ± 0.18 |
| S | FD+H10 | 1.34333 | 0.11015 | 1.34 ± 0.11 |
| S | H14+H10 | 1.64667 | 0.03512 | 1.65 ± 0.04 |
| S | H152+H10 | 1.58667 | 0.09713 | 1.59 ± 0.1 |
| S | H228+H10 | 1.55333 | 0.02309 | 1.55 ± 0.02 |
| TRF | FD+BRS_PF | 1.93000 | 0.05292 | 1.93 ± 0.05 |
| TRF | FD+H10 | 1.91667 | 0.05508 | 1.92 ± 0.06 |
| TRF | H14+H10 | 1.69000 | 0.08185 | 1.69 ± 0.08 |
| TRF | H152+H10 | 1.63667 | 0.08505 | 1.64 ± 0.09 |
| TRF | H228+H10 | 1.66000 | 0.11533 | 1.66 ± 0.12 |
| TRV | FD+BRS_PF | 1.62667 | 0.04509 | 1.63 ± 0.05 |
| TRV | FD+H10 | 1.86000 | 0.02000 | 1.86 ± 0.02 |
| TRV | H14+H10 | 1.39667 | 0.18148 | 1.4 ± 0.18 |
| TRV | H152+H10 | 1.48000 | 0.06083 | 1.48 ± 0.06 |
| TRV | H228+H10 | 1.60000 | 0.17349 | 1.6 ± 0.17 |
VitC
results <- summ_qf(dframe_qf_ufscar,variable_name = "VitC")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| RegHidrico | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| S | FD+BRS_PF | 28.54000 | 0.60622 | 28.54 ± 0.61 |
| S | FD+H10 | 35.68000 | 1.14582 | 35.68 ± 1.15 |
| S | H14+H10 | 32.45667 | 3.28649 | 32.46 ± 3.29 |
| S | H152+H10 | 37.78333 | 1.63102 | 37.78 ± 1.63 |
| S | H228+H10 | 33.36000 | 1.15763 | 33.36 ± 1.16 |
| TRF | FD+BRS_PF | 33.21000 | 1.82371 | 33.21 ± 1.82 |
| TRF | FD+H10 | 32.76000 | 1.84984 | 32.76 ± 1.85 |
| TRF | H14+H10 | 32.50667 | 1.82870 | 32.51 ± 1.83 |
| TRF | H152+H10 | 34.24333 | 0.89131 | 34.24 ± 0.89 |
| TRF | H228+H10 | 34.38333 | 1.21187 | 34.38 ± 1.21 |
| TRV | FD+BRS_PF | 33.45333 | 1.59080 | 33.45 ± 1.59 |
| TRV | FD+H10 | 35.30333 | 0.89891 | 35.3 ± 0.9 |
| TRV | H14+H10 | 38.17667 | 1.34500 | 38.18 ± 1.35 |
| TRV | H152+H10 | 37.89000 | 1.84789 | 37.89 ± 1.85 |
| TRV | H228+H10 | 35.40333 | 2.43097 | 35.4 ± 2.43 |
CCSM
Massa
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Massa")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| Manejo | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| Convencional | BRF_PF+FD | 94.33333 | 11.50362 | 94.33 ± 11.5 |
| Convencional | BRF_PF+H152 | 89.33333 | 10.01665 | 89.33 ± 10.02 |
| Convencional | IAC10+FD | 90.33333 | 4.72582 | 90.33 ± 4.73 |
| Convencional | IAC10+H152 | 86.66667 | 5.13160 | 86.67 ± 5.13 |
| Convencional | IAC5+FD | 88.00000 | 12.12436 | 88 ± 12.12 |
| Convencional | IAC5+H152 | 81.66667 | 2.30940 | 81.67 ± 2.31 |
| Ecologico | BRF_PF+FD | 77.33333 | 3.05505 | 77.33 ± 3.06 |
| Ecologico | BRF_PF+H152 | 83.66667 | 0.57735 | 83.67 ± 0.58 |
| Ecologico | IAC10+FD | 95.00000 | 11.13553 | 95 ± 11.14 |
| Ecologico | IAC10+H152 | 90.33333 | 2.51661 | 90.33 ± 2.52 |
| Ecologico | IAC5+FD | 80.33333 | 3.21455 | 80.33 ± 3.21 |
| Ecologico | IAC5+H152 | 88.00000 | 3.46410 | 88 ± 3.46 |
Altura
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Altura")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| Manejo | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| Convencional | BRF_PF+FD | 6.10000 | 0.40000 | 6.1 ± 0.4 |
| Convencional | BRF_PF+H152 | 5.96667 | 0.25166 | 5.97 ± 0.25 |
| Convencional | IAC10+FD | 6.06667 | 0.15275 | 6.07 ± 0.15 |
| Convencional | IAC10+H152 | 5.56667 | 0.11547 | 5.57 ± 0.12 |
| Convencional | IAC5+FD | 5.90000 | 0.36056 | 5.9 ± 0.36 |
| Convencional | IAC5+H152 | 5.80000 | 0.10000 | 5.8 ± 0.1 |
| Ecologico | BRF_PF+FD | 5.70000 | 0.10000 | 5.7 ± 0.1 |
| Ecologico | BRF_PF+H152 | 5.83333 | 0.05774 | 5.83 ± 0.06 |
| Ecologico | IAC10+FD | 6.13333 | 0.15275 | 6.13 ± 0.15 |
| Ecologico | IAC10+H152 | 6.03333 | 0.11547 | 6.03 ± 0.12 |
| Ecologico | IAC5+FD | 5.90000 | 0.10000 | 5.9 ± 0.1 |
| Ecologico | IAC5+H152 | 6.03333 | 0.05774 | 6.03 ± 0.06 |
Largura
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Largura")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| Manejo | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| Convencional | BRF_PF+FD | 5.60000 | 0.36056 | 5.6 ± 0.36 |
| Convencional | BRF_PF+H152 | 5.40000 | 0.17321 | 5.4 ± 0.17 |
| Convencional | IAC10+FD | 5.40000 | 0.10000 | 5.4 ± 0.1 |
| Convencional | IAC10+H152 | 5.10000 | 0.10000 | 5.1 ± 0.1 |
| Convencional | IAC5+FD | 5.43333 | 0.25166 | 5.43 ± 0.25 |
| Convencional | IAC5+H152 | 5.26667 | 0.05774 | 5.27 ± 0.06 |
| Ecologico | BRF_PF+FD | 5.16667 | 0.11547 | 5.17 ± 0.12 |
| Ecologico | BRF_PF+H152 | 5.26667 | 0.05774 | 5.27 ± 0.06 |
| Ecologico | IAC10+FD | 5.56667 | 0.20817 | 5.57 ± 0.21 |
| Ecologico | IAC10+H152 | 5.23333 | 0.25166 | 5.23 ± 0.25 |
| Ecologico | IAC5+FD | 5.30000 | 0.10000 | 5.3 ± 0.1 |
| Ecologico | IAC5+H152 | 5.43333 | 0.05774 | 5.43 ± 0.06 |
AL
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "AL")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| Manejo | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| Convencional | BRF_PF+FD | 1.09000 | 0.01732 | 1.09 ± 0.02 |
| Convencional | BRF_PF+H152 | 1.10667 | 0.02082 | 1.11 ± 0.02 |
| Convencional | IAC10+FD | 1.12333 | 0.01155 | 1.12 ± 0.01 |
| Convencional | IAC10+H152 | 1.09333 | 0.01155 | 1.09 ± 0.01 |
| Convencional | IAC5+FD | 1.08667 | 0.02082 | 1.09 ± 0.02 |
| Convencional | IAC5+H152 | 1.10000 | 0.01000 | 1.1 ± 0.01 |
| Ecologico | BRF_PF+FD | 1.10667 | 0.03055 | 1.11 ± 0.03 |
| Ecologico | BRF_PF+H152 | 1.10667 | 0.01528 | 1.11 ± 0.02 |
| Ecologico | IAC10+FD | 1.10333 | 0.01155 | 1.1 ± 0.01 |
| Ecologico | IAC10+H152 | 1.15333 | 0.03786 | 1.15 ± 0.04 |
| Ecologico | IAC5+FD | 1.11333 | 0.02082 | 1.11 ± 0.02 |
| Ecologico | IAC5+H152 | 1.11000 | 0.02000 | 1.11 ± 0.02 |
Suco
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Suco")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| Manejo | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| Convencional | BRF_PF+FD | 50.00000 | 2.66646 | 50 ± 2.67 |
| Convencional | BRF_PF+H152 | 48.80000 | 1.57162 | 48.8 ± 1.57 |
| Convencional | IAC10+FD | 48.30000 | 0.55678 | 48.3 ± 0.56 |
| Convencional | IAC10+H152 | 47.23333 | 2.61024 | 47.23 ± 2.61 |
| Convencional | IAC5+FD | 48.76667 | 1.51438 | 48.77 ± 1.51 |
| Convencional | IAC5+H152 | 50.63333 | 1.25033 | 50.63 ± 1.25 |
| Ecologico | BRF_PF+FD | 51.66667 | 2.72274 | 51.67 ± 2.72 |
| Ecologico | BRF_PF+H152 | 51.00000 | 0.34641 | 51 ± 0.35 |
| Ecologico | IAC10+FD | 49.13333 | 1.68028 | 49.13 ± 1.68 |
| Ecologico | IAC10+H152 | 45.43333 | 3.35609 | 45.43 ± 3.36 |
| Ecologico | IAC5+FD | 51.03333 | 1.51767 | 51.03 ± 1.52 |
| Ecologico | IAC5+H152 | 49.20000 | 1.76918 | 49.2 ± 1.77 |
Acidez
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Acidez")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| Manejo | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| Convencional | BRF_PF+FD | 5.79277 | 0.21457 | 5.79 ± 0.21 |
| Convencional | BRF_PF+H152 | 4.73127 | 1.12089 | 4.73 ± 1.12 |
| Convencional | IAC10+FD | 5.78470 | 0.14056 | 5.78 ± 0.14 |
| Convencional | IAC10+H152 | 4.64587 | 0.14255 | 4.65 ± 0.14 |
| Convencional | IAC5+FD | 5.50670 | 0.55673 | 5.51 ± 0.56 |
| Convencional | IAC5+H152 | 5.28863 | 0.23447 | 5.29 ± 0.23 |
| Ecologico | BRF_PF+FD | 5.43200 | 0.25991 | 5.43 ± 0.26 |
| Ecologico | BRF_PF+H152 | 5.71203 | 0.41669 | 5.71 ± 0.42 |
| Ecologico | IAC10+FD | 5.43063 | 0.73447 | 5.43 ± 0.73 |
| Ecologico | IAC10+H152 | 5.57807 | 0.10602 | 5.58 ± 0.11 |
| Ecologico | IAC5+FD | 5.59083 | 0.15735 | 5.59 ± 0.16 |
| Ecologico | IAC5+H152 | 5.64403 | 0.08699 | 5.64 ± 0.09 |
BrixCorr
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "BrixCorr")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plo ## NULLTabela Resumo
results$dframe_summ %>% kable_cust()| Manejo | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| Convencional | BRF_PF+FD | 9.03333 | 0.15275 | 9.03 ± 0.15 |
| Convencional | BRF_PF+H152 | 9.03333 | 0.15275 | 9.03 ± 0.15 |
| Convencional | IAC10+FD | 9.06667 | 0.15275 | 9.07 ± 0.15 |
| Convencional | IAC10+H152 | 9.00000 | 0.30000 | 9 ± 0.3 |
| Convencional | IAC5+FD | 9.26667 | 0.05774 | 9.27 ± 0.06 |
| Convencional | IAC5+H152 | 9.30000 | 0.17321 | 9.3 ± 0.17 |
| Ecologico | BRF_PF+FD | 9.26667 | 0.15275 | 9.27 ± 0.15 |
| Ecologico | BRF_PF+H152 | 9.33333 | 0.15275 | 9.33 ± 0.15 |
| Ecologico | IAC10+FD | 9.13333 | 0.15275 | 9.13 ± 0.15 |
| Ecologico | IAC10+H152 | 8.86667 | 0.05774 | 8.87 ± 0.06 |
| Ecologico | IAC5+FD | 9.40000 | 0.10000 | 9.4 ± 0.1 |
| Ecologico | IAC5+H152 | 8.96667 | 0.15275 | 8.97 ± 0.15 |
Ratio
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "Ratio")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| Manejo | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| Convencional | BRF_PF+FD | 1.56667 | 0.05774 | 1.57 ± 0.06 |
| Convencional | BRF_PF+H152 | 1.96667 | 0.45092 | 1.97 ± 0.45 |
| Convencional | IAC10+FD | 1.56667 | 0.05774 | 1.57 ± 0.06 |
| Convencional | IAC10+H152 | 1.93333 | 0.05774 | 1.93 ± 0.06 |
| Convencional | IAC5+FD | 1.66667 | 0.15275 | 1.67 ± 0.15 |
| Convencional | IAC5+H152 | 1.76667 | 0.11547 | 1.77 ± 0.12 |
| Ecologico | BRF_PF+FD | 1.73333 | 0.05774 | 1.73 ± 0.06 |
| Ecologico | BRF_PF+H152 | 1.63333 | 0.11547 | 1.63 ± 0.12 |
| Ecologico | IAC10+FD | 1.70000 | 0.26458 | 1.7 ± 0.26 |
| Ecologico | IAC10+H152 | 1.56667 | 0.05774 | 1.57 ± 0.06 |
| Ecologico | IAC5+FD | 1.66667 | 0.05774 | 1.67 ± 0.06 |
| Ecologico | IAC5+H152 | 1.56667 | 0.05774 | 1.57 ± 0.06 |
IT
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "IT")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| Manejo | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| Convencional | BRF_PF+FD | 1.84667 | 0.06807 | 1.85 ± 0.07 |
| Convencional | BRF_PF+H152 | 1.80333 | 0.03215 | 1.8 ± 0.03 |
| Convencional | IAC10+FD | 1.79000 | 0.02646 | 1.79 ± 0.03 |
| Convencional | IAC10+H152 | 1.73333 | 0.04726 | 1.73 ± 0.05 |
| Convencional | IAC5+FD | 1.84000 | 0.05292 | 1.84 ± 0.05 |
| Convencional | IAC5+H152 | 1.92000 | 0.03606 | 1.92 ± 0.04 |
| Ecologico | BRF_PF+FD | 1.94667 | 0.07506 | 1.95 ± 0.08 |
| Ecologico | BRF_PF+H152 | 1.94333 | 0.01528 | 1.94 ± 0.02 |
| Ecologico | IAC10+FD | 1.83667 | 0.03215 | 1.84 ± 0.03 |
| Ecologico | IAC10+H152 | 1.64000 | 0.12124 | 1.64 ± 0.12 |
| Ecologico | IAC5+FD | 1.95333 | 0.04163 | 1.95 ± 0.04 |
| Ecologico | IAC5+H152 | 1.79667 | 0.09018 | 1.8 ± 0.09 |
VitC
results <- summ_qf_ccsm(dframe_qf_ccsm,variable_name = "VitC")Gráfico de Barras
results$barp
results$barp_2
Gráfico De Barras com erro
results$barp_error
results$barp_error_2
Gráfico de Pontos + Média e Erro
results$mean_desv_plot
results$mean_desv_plot_2
Tabela Resumo
results$dframe_summ %>% kable_cust()| Manejo | Tratamento | Mean | SD | Mean ± SD |
|---|---|---|---|---|
| Convencional | BRF_PF+FD | 34.15667 | 2.19805 | 34.16 ± 2.2 |
| Convencional | BRF_PF+H152 | 37.52333 | 1.58418 | 37.52 ± 1.58 |
| Convencional | IAC10+FD | 34.25333 | 1.33121 | 34.25 ± 1.33 |
| Convencional | IAC10+H152 | 33.71000 | 2.57874 | 33.71 ± 2.58 |
| Convencional | IAC5+FD | 32.83667 | 2.87114 | 32.84 ± 2.87 |
| Convencional | IAC5+H152 | 33.98333 | 1.23678 | 33.98 ± 1.24 |
| Ecologico | BRF_PF+FD | 31.96000 | 1.10585 | 31.96 ± 1.11 |
| Ecologico | BRF_PF+H152 | 33.32333 | 3.61777 | 33.32 ± 3.62 |
| Ecologico | IAC10+FD | 34.23000 | 2.17690 | 34.23 ± 2.18 |
| Ecologico | IAC10+H152 | 36.66667 | 0.76846 | 36.67 ± 0.77 |
| Ecologico | IAC5+FD | 31.53667 | 1.24340 | 31.54 ± 1.24 |
| Ecologico | IAC5+H152 | 33.23333 | 1.44901 | 33.23 ± 1.45 |